Numerical Computation of Thermoelectric and Thermomagnetic Effects 
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Abstract 

Phenomenological equations describing the Seebeck, Hall, 
.Nernst, Peltier, Ettingshausen, and Righi-Leduc effects are 
■numerically solved for the temperature, electric current, and 
electrochemical potential distributions of semiconductors un- 
der magnetic field. The results are compared to experiments. 

1 Introduction 

It is well known that geometry of samples (e.g., ratio of 
length to width for rectangular samples) affects magnetoresis- 
tance. Short samples become more electrically resistant under 
magnetic field than long ones (see Fig. ^below). 

1 For the Seebeck and the Nernst effects, experimental evi- 
■dences for such geometric contribution are not so clear. Ertl 
measured Bi-Sb alloy samples of various lengths, and showed 
that longer samples exhibit greater Seebeck coefficients. Mea- 
surement of Seebeck coefficient by Ikeda et al however, 
is not easy to summarize, but their Nernst coefficient of wider 
X"fat bridge") sample (Fig. under 4 Tesla was about 12 per- 
cent smaller than that of narrower one (Fig. 

Though analytic solutions exist for a limited class of the Hall 
effect [|||, numerical computation is necessary to explain these 
results in general. We developed a two-dimensional numeri- 
cal simulation code based on phenomenological equations gov- 
|erning thermoelectric and thermomagnetic effects. Section ^ 
■summarizes the basic equations. Sections |^-|^ describe the nu- 
,merical algorithm. Sections ^ and ^ summarizes the results and 
■discuss consequences. 

2 Phenomenological Equations 

The phenomemological equations governing_ thermoelectric 
and thermomagnetic phenomena of solids are | 
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-V(i) = pJ + aVT + RB X J + NB xVT (1) 
q = (t)J - kVT + aTJ + NTB x J + kMB x VT (2) 

where (j) is the electrochemical potential per unit charge,[] p the 
(isothermal) electric resistivity, J the electric current density, 
a the (isothermal) Seebeck coefficient, T the temperature, R 
the (isothermal) Hall coefficient, B the magnetic flux density, 
N the (isothermal transverse) Nernst coefficient, q the energy 
flux density, k the (isothermal) thermal conductivity, and AI the 

'For electrons, 4' = 0clcctrosiaiic + /^n/(— e), where is the chemical 
potential of electrons with charge — e < 0. Similarly, cf> = </<clcctiostatic + 
fip/ (+e) for holes with charge +e > 0. If the system is close to equilibrium, 
the (/)'s for each species of carriers differ very little, so we shall not distinguish 
the <A's and the J's for different carriers. 



Figure 1 : Example of discretization. This figure corresponds to 
the "bridge" shape of Ikeda et al [^]. 

Righi-Leduc coefficient.^ In Eq. ([l|), the last three terms repre- 
sents the Seebeck, the Hall, and the Nernst effects, respectively. 
In Eq. (||), the last three terms are responsible for the Peltier/ 
Thomson, the Ettingshausen, and the Righi-Leduc effects, re- 
spectively. 

In what follows we assume that (1) the system is in steady 
state, so that V ■ J = V ■ q = holds; (2) the external mag- 
netic field B is independent of the position,^ and is along the 
z-direction; (3) the electric current J has no z-component; (4) 
the temperature T is independent of the z-coordinate; (5) the 
conductor is homogenious, so that all of the transport coeffi- 
cients (a, p, K, R, N, AI) are functions of temperature T alone. 

3 Overview of the Method 

Our aim is to calculate T, (p, and J distributions for rectangu- 
lar and irregular-shaped semiconductor samples such as shown 
in Fig. |l[ We discretize position by constructing a grid with 
square meshes of size h x h. We consider T and (j> on each 
grid points (corners of the meshes), and J on each side of the 
meshes, as shown in Fig. ^. 

After setting up suitable initial values for T, cp, and J, we 
proceed as follows: 

1. On each grid point update Tij from the discretized 
Poisson equation for temperature (see Section ^ below), 
assuming all the other quantities fixed. 

2. On each line segment connecting adjacent grid points, up- 
date Jij (see Section ^ below), assuming all the other 
quantities fixed. 



^Note that the "Leduc-Righi coefficient" L of Landau, Lifshitz, and 
Pitaevskii |^ corresponds to our kM. 

''Rigorously speaking, external field B gives rise to electric current J, 
which in turn modifies the field B according to the Maxwell equation V X B = 
fij + ^/idE/dt. (In steady .state dE/dt = 0.) In practice, however, magnetic 
permeability ^ fio = Att X 10^^ is so small that V X B can be safely 
neglected. 
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Figure 2: Grid point and its four adjacent points. Temper- 
ature Tij and potential 0^ are specified on point whereas 
electric current J^j is specified along the line segment connect- 
ing two points and {i + l,j). 



The normal derivative, say VyT, can be derived from Eq. (||) 
if there is no energy and current transfer through the boundary 

(Qy ^ Jy^ 0), 



(7) 



5 Electric Current Calculation 

Given T and (f), the electric current density J can be com- 
puted from Eq. (n|), which can be written as 



pJ + RB X J = -V4> 
or, in coordinate components 



aVT - NB X VT 



(8) 
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Hence 



RB, 
P 



NB.VyT 



-V. 



-aVyT-NB.V^T 



(9) 



3. On each grid point, update (pij so as to satisfy the continu- 
ity equation for at the point (see Section ^ below). 

4. Go to Step 1. 

4 Temperature Updates 

The Poisson equation for temperature can be derived by tak- 
ing the divergence of Eq. (0) and using Eq. (|l]): 



y2y 



P 



NTB"^ dR 
~p If 



(VT) • J 



fNTB^ dN 
\ P 



dT 



(3) 



This second-order partial differential equation for T, which 
we shall abbreviate as V^r(x, y) = F(x,y), can be discretized 
as follows: If the point is not on the boundary. 



rjincw 
ij 



-F (4) 



(with suitable modification to accelerate convergence). On the 
boundary, either is given (Dirichlet conditions), or deriva- 
tives of T in the direction normal to the boundary, V„T(a;, y), 
is given (Neumann conditions). In the latter case, if the grid 
point (i, j) is on the boundary which is along the x-direction 
such that point [i — is outside of the sample, the normal 
derivative VyT = dT/dy should be given, and the update for- 
mula is 



Tnncw 

^ij — 



T, 



T, 



_^ —-^F (5) 



Similarly, at a corner point such that points {i — and {i, j — 
1) are outside of the boundary, the update formula becomes 

= A -y^ 



J. 
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RB, 
P 

aVyT - NB,V,T 



NB.VyT 



(10) 



On the boundary, if Jy = holds, then 

Jcc^i~V.q}-aV.,T + NB,yyT)/p (11) 
Furthermore, if Qy — also holds, substitution of (0) into ( pT| ) 



gives 



p 



-V^cj) -{a- NMBDV^T - 



(12) 



Solving for J^, we have 

Jx = (-V,(/) - (a - NMBDV^T) 



N^B^,T 



(13) 



6 Potential Updates 

Instead of using a lengthy Poisson equation for the electro- 
chemical potential (j) that can be derived from Eq. ([l|), we base 
our calculation of ((> on the continuity equation of electric cur- 
rent, V • J = 0. 

We start with Eq. (p^ which has the form 



1 



R'^B 



2 {-P'^x<t> - RB,\JyC, 



terms independent of (14) 



and similarly for Jy. As was discussed earlier, the discretized 
value Jf^ is taken along the line segment connecting two points 



with potentials (f>ij and ( 



Along this line segment, we can 



approximate potential derivatives by Vj 
and Vyd 



^i,,_i)/(4/j).Af- 
ter these derivatives are substituted, the equation for Jf^ above 
has the 0^ dependence: 



h{p^+R^B2) 



terms independent of < 



(15) 




0.5 1 

Magnetic Induction (Tesia) 

Figure 3: Simulation results of magnetoresistance for intrin- 
sic InSb semiconductor with different geometry: length (x- 
direction) x width (y-direction) — 10 mm x 1mm, 10 mm 
X 10 mm, and 4 mm x 12 mm, with negligible thickness (z- 
direction). External magnetic field is along the z-direction. 



Note that J[j is the current outgoing from the point in the 
positive x-direction. As can be seen from Fig. ^ there are four 



such expressions outgoing from point (i, j): Jf,, Jf-, 



-Ij' 

When these four expressions are added, we arrive at 



J° 



h{p^+R'^B^) 
Now, if we replace the value of 



terms independent of ( 



by 



/now 



j9-^hip' + R'B')/{Ap) 



(16) 



(17) 



the right-hand side of Eq. (16) will vanish. We use Eq. ( |17[ ) and 
similar ones with suitable modification for grid points on the 
boundary, to update 0y . 

7 Results 

We conducted numerical computations on intrinsic indium 
antimonide (InSb) semiconductor samples with the following 
properties near 300 K. 



p = 


8 X IQS 


Q m 


R = 


(_5.6e-0-034(T-273) _ g.Q) x IQ-^ 


m3 A^^ s^^ 


a = 


(-3.2 + 0.01(r - 273)) X 10-"^ 




N = 


(-5.7e-(^-273)/65 _ 3.2) X IQ-^ 




M = 


5 X 10-2 




K = 


1.4 X 105T-1-65 


WK-im-i 



(18) 

These values are not meant to be good fits to measurements. 
They are only shown here as example inputs to our code. (In 
fact, p, a, and N are rough fits to measurements near 300 K by 
Ikeda and others [|[|[ ^, but the conditions are not uniform: N 
is measured under 4 Tesla, whereas a under no magnetic field.) 

On the basis of these values, we computed magnetoresistance 
and Seebeck and Nernst effects with various sample geometry. 

The magnetoresistance results (Fig. |]) are in good agreement 
with experiment Fig. 4.8 of Seeger [^). 

Seebeck and Nernst coefficients should vary much with the 
magnetic field, so the results for these coefficients are to be 
compared with those of different geometry under the same mag- 
netic field. 
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Figure 4: Simulation results of apparent Seebeck coefficient 
for intrinsic InSb samples with different geometry. Apparent 
Seebeck coefficient is defined to be acfr = Ax(j>/ A^T, where 
Ax(f> is longitudinal potential difference, and A^T is longitu- 
dinal temperature difference. Solid lines: Current leads are of 
neghgible widths. Dashed lines: Current leads are as wide as 
the sample widths. 



As Fig. ^ shows, the Seebeck effect is not sensitive to ge- 
ometry if current leads that measure longitudinal voltage dif- 
ference are narrow enough. On the other hand, if current leads 
are as wide as the sample width, Seebeck effect degrades and 
can even change sign for short samples. This tendency ex- 
plains some experimental evidence for the size dependence of 
magneto-Seebeck effect [|l|. 

Similar tendency can be seen in Fig. ^ for the Nernst ef- 
fect. In this case, however, geometry effect exists even with 
narrow current leads. This is because the Righi-Leduc effect 
causes transverse (y-direction) temperature gradient that is pro- 
portional to the magnetic field B. This transverse temperature 
gradient in turn causes transverse voltage gradient by the See- 
beck effect. When divided by B, this transverse voltage gives a 
nearly constant bias to the apparent Nernst coefficient. 

This effect partly explains what Ikeda et al [|[] found out: Un- 
der 4 Tesla the apparent Nernst coefficient for their "fat bridge" 
(wide with legs) sample (Fig. 7 ) is about 12 percent smaller than 
those of narrower one (Fig. 8), whereas our simulation gives 
only 7 to 8 percent smaller coefficient. Though they are careful 
to make their current leads narrow, inevitable finite widths of 
the leads might explain further effect. 

8 Discussion 

As can be seen from the contour maps, the potential distribu- 
tions of thermomagnetic samples are rather complicated. They 
change dramatically with different sample geometry. More- 
over, if we attach current leads of finite widths to the cold and 
hot edges, much of the transverse voltage gradient is shorted 
out, resulting in quite a different potential distribution. To cor- 
rect for such a bias, careful numerical calculation is necessary. 

Inspection of Fig. ^ (upper) and Fig. ^ shows that the best 
way to generate electricity from wide samples under magnetic 
field is to attach current leads to the bottom left and the top 
right corners. In this way, we can utilize both the Seebeck and 
the Nernst effects. 

Figs. ^ and|^ also show that at the far end of a small "leg," 
potential gradient along the boundary decreases by an order of 
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Figure 5: Simulation results of apparent Nernst coefficient 
for intrinsic InSb samples with different geometry. Apparent 
Nernst coefficient is defined to be A^cff = LAy(j)/BWAxT, 
where Ayip is transverse potential difference, AxT longitudi- 
nal temperature difference, L length (in x-direction), W width 
(in y-direction), and B magnetic induction. Solid lines: Cur- 
rent leads are of negligible widths. Dashed lines: Current leads 
are as wide as the sample widths. 




Figure 6: Contour maps of calculated potential distribution for 
a 12 mm x 4 mm sample in 4-Tesla magnetic field. Upper: Cur- 
rent leads are of negligible width. Lower: Current leads are as 
wide as the sample width. The left and the right edges are kept 
at 273 K and 373 K, respectively. 
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Figure 8: Contour map of calculated potential distribution for 
"bridge" sample of Ikeda et al. 



magnitude or more, thus making measurement of Nernst coef- 
ficients less prone to the widths of current leads. 
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Figure 7: Contour map of calculated potential distribution for 
"fat bridge" sample of Ikeda et al [^. 



